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ABSTRACT 

We present Keck spectroscopic measurements of the millisecond pulsar binary J2215+5135. These data 
indicate a neutron-star (NS) mass Mns = 1.6 Mq, much less than previously estimated. The pulsar heats 
the companion face to Tjj ps 9000 K; the large heating efficiency may be mediated by the intrabinary shock 
dominating the X-ray light curve. At the best-fit inclination i = 88.8°, the pulsar should be eclipsed. We find 
weak evidence for such eclipses in the pulsed gamma-rays; an improved radio ephemeris allows use of up to 
5 times more Fermi-LAT gamma-ray photons for a definitive test of this picture. If confirmed, the gamma-ray 
eclipse provides a novel probe of the dense companion wind and the pulsar magnetosphere. 

Subject headings: gamma rays: stars — pulsars: general 


1. INTRODUCTION 

PSR J2215+5135 is a recycled millisecond pulsar (MSP) 
with P = 2.6 ms and E = 7 . 4/45 x 10 34 ergs _1 (where 
/45 is the neutron star [NS] moment of inertia in units of 
10 45 gem 1 2 ), in a P 5 = 4.14hr orbit with a heated low-mass 
stellar companion. This “redback” (RB) system was dis¬ 
covered in 350 MHz Green Bank Telescope (GBT) observa¬ 
tions (reported by Ray et al. 2012) of an unidentified Fermi- 
LAT (Large Area Telescope. I Atwood et al.1:2009 ) gamma-ray 
source. The radio dispersion measure (69.2 cm ~ 3 pc) pro¬ 
vides a distance estimate 3 kpc, and this relatively close 
and bright system can be subject to detailed study. 

Schroeder & Halpern (2014, hereafter SH14) obtained 
high-quality BVR light curves of the companion over many 
orbits. It varies from V ~ 18.7 to 20.2 mag and thus is well 
measured with modest-aper ture telescopes. Fitting thes e light 
curves with the ELC code (lOrosz & Hausch i Idt]120001) and a 
photometr y tab le generated from the PHOENIX model atmo¬ 
spheres dHusser et al.ll2013 ). SHI4 inferred a low system in¬ 
clination and large pulsar mass, M\ = 2A5+qH Mq. This 
is potentially very important, since it would be one of the 
highest-known NS masses. Such objects can provide strong 
constraints on the equation of state at supern uclear densities 
jLattimer & Prakashll2007b ISteiner et al.ll2013T> . 

However, their fit showed several peculiarities. The un¬ 
heated “night” face of the companion was found to have 
T c = 3790^25 K and the heated “day” side T D = 4899K 
in the ELC fits, in poor agreement with the observed bluer 
colors. The light curve seemed shifted by A cf> « —0.01 with 
respect to the radio-pulse ephemeris. SH14 suggest that the 
blue colors could be a hot quiescent disk, while the phase shift 
may be an effect of an intrabinary shock. Also, iGentile et al.l 
(120141) observed the system in the X-rays with CXO, finding 
an X-ray minimum near orbital phase <f> « 0.25 (pulsar su¬ 
perior conjunction, radio eclipse), which they interpret as due 
to variable obscuration of emission from an intrabinary shock 
around the companion. 

To probe these peculiarities, we have obtained Keck spectra 
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across the binary orbit. Our data show no evidence for disk 
or shock emission lines, but imply much higher system in¬ 
clination, with higher companion temperatures and stronger 
pulsar heating. This lowers the inferred pulsar mass, but 
raises the possibility that the pulsed magnetospheric emission 
is eclipsed by the companion star. Gamma-ray data currently 
show limited evidence for such eclipse, but additional obser¬ 
vations can test this hypothesis. 

2. KECK SPECTROSCOPY 

The system was observed on three occasions. In each case 
we aligned the long slit at position angle 76° (N through E) 
to include two bracketing stars for monitoring purposes (see 
the finder chart in Breton et al. 2013). The R = 17.6 mag 
star at USNO-B1 J2000 position a = 22 hr 15 m 32.031 s , 5 = 
+51°35'35.46" (Cl) appears to show radial-velocity (RV) 
variations between epochs, so we use the R = 18.1 mag 
star 8" away at J2000 position a = 22 hr 15 m 33.498 s , 5 = 
+51°35 , 38.88" (C2) as our reference object. 

First, using the DEep Imaging Multi-Object Spectrograph 
(DEIMOS; Faber et al. 2003) on Keck-II, we obtained 
7 x 600 s exposures covering 4450-9060 A on 2014 October 
2 (UT dates are used throughout this paper; MJD 56932) with 
~ 4.7 A resolution through a l"-wide slit. These data cov¬ 
ered 4>b = 0.685-1.025, where <pB = 0.75 is optical maxi¬ 
mum (“Noon”) and the pulsar ascending node is at <f>B = 0. 
The radio ephemeris defines Tasc = 55186.164485831MJD 
and Pb = 0.172502104907d. For improved blue cover¬ 
age, we observed again on 2014 October 24 (MJD 56954) by 
taking 6 x 600 s exposures with the Low Resolution Imag¬ 
ing Spectrometer (LRIS; Oke et al. 1995) on Keck-I us¬ 
ing the 1" slit and the 5600 A dichroic splitter. The blue- 
camera 600/4000 grism provided ~ 4 A resolution, while 
the red-camera 400/8500 grating delivered ~ 7 A resolution. 
These data covered wavelengths ~ 3200-10.260A and or¬ 
bital phases 0.501 < 4>b < 0.723. Finally, we observed with 
the same setup on 2014 November 20 (MJD 56981), covering 
orbital phase 0.202 < 4>b < 0.381 by taking 5 x 600 s ex¬ 
posures. Calibration spectra of blue (BD+28°4211) and red 
(BD+17°4708) spectrophotometric standards were obtained 
during each run. The spectra were subject to standard IRAF 
reductions, including optimal extraction. 

Figure 1 shows the reduced, calibrated LRIS spectra at three 
orbital phases. The companion is significantly hotter than 
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Fig. 1.— Keck-I/LRIS spectra of J2215+5135, showing maximum bright¬ 
ness (“Day”; <f>g = 0.73), “Dawn” (o b = 0.50), and “Night” (<Pb = 
0.20). Each spectrum is a single 600 s exposure. Broad features near 5500 A 
are artifacts from incomplete correction for the dichroic response. 

the SH14 ELC solution, with an effective spectral class A2 V 
near optical maximum brightness (T e ff ss 8970 K) and GOV 
near minimum (T e g ss 6030 K), as determined by line-ratio 
comparison with spectral standards. No strong emission lines 
or nonthermal component are seen, although a weak variable 
blue continuum peaking at L/> « 0.75 is not at present ex¬ 
cluded. 

The spectra show a wealth of stellar absorption features, 
and so we measure the RV variation using the IRAF XCASO 
script dKurtz & Minkl 119981) . The templates were drawn 
from the Indo-US library of Coude Feed spectra (Vald es et al.1 
120041) . selecting dwarf stars with known RV from a range of 
well-determined spectral classes. The cross-correlation statis¬ 
tic was R « 30-60 for the day-phase October 24 LRIS data 
and R « 20-25 for the night-phase November 20 spectra. 
With the more-limited blue coverage on October 2, the cross 
correlations were weaker, giving R « 10-20 and larger RV 
uncertainties. Since Cl shows ~ 50kms -1 RV variation 
between epochs, we relied on C2 for our RV tie. This star 
matched best to a K0 V spectrum. Cross correlation delivered 
velocity uncertainties of ~ 2.5kms _1 (LRIS) and ~ 3.5— 
4.0 kms -1 (DEIMOS). To remove residual drift in the wave¬ 
length solution, we corrected the pulsar companion veloci¬ 
ties to bring the measured C2 velocities into agreement. The 
C2 measurement errors, plus systematic 3 km s -1 errors esti¬ 
mated from the variance between cross correlations with dif¬ 
ferent templates, were added in quadrature to the individual 
RV uncertainties to produce the final RV uncertainty. These 
RVs and errors are plotted in Figure [2] as a function of orbital 
phase computed from the barycentered times of the exposure 
midpoints. Lacking velocity-standard observations spanning 
the companion spectral class during each run, we may have 
additional systematic uncertainty in the absolute RV. 

3. LIGHT-CURVE COMPARISON AND MODELING CONSTRAINTS 

To compare with the light-curve data, we plot the phased 
BVR magnitudes and uncertainties published in Figure 2 of 
SHI4; we ignored upper-limit points (generally from epochs 



Fig. 2.— Top: BV R MDM photometry l lSchroeder & Halpernll2014l) . Bot- 
tom: Keck DEIMOS/LRIS radial velocities. Two periods are shown, with 
data and the large-z (“HiT”) ELC model in the first cycle. The second period 
shows the model residuals, multiplied by a factor of 3. In the velocity panel, 
we also show the residuals after a best-fit sinusoid is removed from both the 
model and the data. 

of poor photometry), extracting 103 B , 55 V, and 113 7? mag¬ 
nitudes. These data and our new RV measurements were an¬ 
alyzed with two popular light-curve modeling programs: the 
ELC code and the ICARUS code (Breton et al. 2011). For 
both we can run in “MSP” mode, with the projected pul¬ 
sar orbit x\ = aisini = 0.4681414331t-s from the radio 
ephemeris. Both codes model pulsar heating of the compan¬ 
ion as illumination by a point (X-ray) source, generating filter- 
specific light-curve models and estimating binary-system pa¬ 
rameters via model fits. The temperature model is determined 
by the underlying temperature of the star (actually T c of the 
unheated “night” face), and a heating flux denoted L //. The 
observed light curves are then sensitive to the orbital inclina¬ 
tion i and the mass ratio q = M^/M c . The heating power 
can be related to the effective temperature of the heated face 
T d by 

L h = (Tp - T^)47r[iri(l + g)] 2 <r/sin 2 * (1) 

(effective albedo = 0). The Roche lobe filling factor f c is also 
relevant; as concluded by SH14, ELC fits require f c « 1. 
It is also large (f c > 0.86) in the ICARUS fits. SH14 find 
that the heating center is phase shifted. We concur, but find 
Acj) = —0.0089, smaller than the A <j> = —0.0140 ± 0.0005 
of their Table 2. 

The two codes have some model differences. For the ELC 
code, one can simultaneously fit the observed RV points, giv¬ 
ing additional sensitivity to q and direct estimates of the com¬ 
ponent masses. However, the code fits normalized light curves 
for each band separately, ignoring the instantaneous colors 
and observed magnitude. When searching high-quality data 
for the minimum near a correct physical model, this normal¬ 
ization should not affect the fit values and insulates the results 
from systematic interband photometry errors. The ICARUS 
code, in contrast, fits the observed magnitudes directly, using 
the instantaneous colors and providing estimates of the dis¬ 
tance modulus and the extinction. 
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TABLE 1 

Binary Model Fits 


Param. 

SH14 

pa 

P/S b 

HiT c 

ICARUS 

i 

51-6±|I 

51.8 

52.0 

88.8 

75-90 

<? 

6.20 ± 0.25 

7.50 

7.35 

6.89 

(0.87-0.99f 

T c ( K) 

3790+25 

3780 

3780 

6220 

6637±12 

Log (L h ) 

33.79 

33.79 

33.77 

34.71 

34.76 

X 2 /v 

[3.29] e 

3.25 

2.98 

4.36 

5.05 

Mns 

z - 4O -0.11 

4.14 

3.88 

1.59 

- 

M c 

0.396 ± 0.045 

0.55 

0.53 

0.23 

— 


a P: ELC photometric fit 

b P/S:ELC photometric/spectral fit (“Best” in Figure 3.) 
c HiT ELC photometric/spectral fit w/ T c > 6000 K constraint. 
d Roche fill factor. f c increases with v; all / c > 0.9 imply i ss 90°. 
e Our fit value; SH14 report \' 2 /v ss 1.5. 

For the ELC modeling, we wish to compare with the SH14 
results, and so we used the same color table (kindly shared by 
J. Schroeder) generated from the PHOENIX atmosphere mod¬ 
els dHusser et al.ll2013l) with an extension to T e g > 10,000 K 
from the ATLAS9 atmospheres (ICastelli & Kuruczl 120041) . 
For the ICARUS code we collected Harris BVR color tables 
from the PHOENIX models tabulated at the Spanish Virtual 
Observatory (svo2.cab.inta-csis.es). Our results are summa¬ 
rized in Table [Hand Figure [3] 

Our ELC fit to the BVR magnitudes finds parameters sim¬ 
ilar to those of SH14, except that their solution lies to the side 
of a broad minimum; the global minimum lies at consider¬ 
ably higher q ss 7.5, with x 2 /u = 3.25 (the minimum in 
their estimated ler contours is x 2 / u = 3.27). These x 2 fv 
values are about twice those reported by SH14, but the model 
light curves and residual plots look quite similar to their re¬ 
sults across this range of q. These ELC photometry fits give 
much lower temperatures and heating powers than allowed by 
our Keck spectra. The ICARUS photometry fit gives a dra¬ 
matically different picture, finding T c « 6600 K, much more 
consistent with the spectral estimate. Lh is also much larger. 
The fit does not strongly constrain q. The best fits prefer large 
extinctions; we fix Ay = 0. 62 mag, the maximum in this di¬ 
rection from the iSchlaflv & Finkbeineii (1201 11) models. The 
distance modulus is 12.98mag or d ~ 3.9kpc, not inconsis¬ 
tent with the dispersion-measure estimate. 

We next explored ELC fits of the combined photometry and 
spectroscopy RV datasets. Since there are only 18 velocity 
points, these could be overwhelmed by the photometric data. 
Accordingly, we fit using the RV points 5 times each, effec¬ 
tively increasing the weights to compete with the photometry 
sets. Testing with only 1 and as many as 10 RV sets revealed 
that the fit results are insensitive to this choice, but the fitting 
converged more rapidly to the minima with the increased RV 
weights. 

Figure 3 illustrates the results from our combined photom- 
etry/RV fits, with two panels showing the projected x 2 from 
a grid of models in four fit parameters. The points are scaled 
so that large sizes correspond to lower x 2 - With a minimum 
X 2 v, we must conclude that the models are inadequate 
(or the error bars are underestimated). Thus, no particular 
fit value or formal uncertainty should be believed. Neverthe¬ 
less, the global picture of a pulsar-heated companion is surely 
correct, and the distribution of best (albeit inadequate) fits in 
parameter space gives a guide to future precise solutions. 

First, we note the strong correlation of the “good” solutions 


in the T c -Lh plane. Recalling that ELC fits the normalized 
light curves, we can understand that this simply preserves the 
height of the light-curve peak; with higher T c , larger L u is 
needed for the same maximum. The best light-curve fits are 
indeed for low T c solutions, as found by SHI4. However there 
is a shallower minimum at 6000 K < T c < 7000 K, con¬ 
sistent with the colors and spectroscopic temperatures. This 
local minimum is not evident when fitting only photometry. 
The RV data require q substantially larger than in the pure 
photometry fits, but i remains poorly constrained, with a wide 
swath of moderate x 2 displaying many local minima. The 
small SH14 q = 6.2 value is excluded, but the global mini¬ 
mum i « 52° is quite similar to the best-fit photometry-only 
solution. This solution has an unphysically large pulsar mass 
-Mns = 3.88 Mg. 

However, if we restrict attention to solutions with T c > 
6000 K, consistent with the spectroscopic classification 
(square dots), these form a minimum at large inclination i, and 
q « 6.9. This shallower but well-defined minimum, shown in 
detail in the upper-panel inset of Figure 3, has a much more 
modest Mns = 1.59 Mg. In Figure 2, we show the light- 
curve and RV models for this solution, and the fit residuals. 
Note the systematic residuals over optical maximum, espe¬ 
cially in V and R. This indicates an incorrect heating model. 
The correct heating pattern (likely mediated by an intrabinary 
shock; see below) should have a surface temperature distribu¬ 
tion allowing a better match to the light-curve peak. The RV 
residuals show similar systematic departures. We also plot the 
RV points and model after the same best-fit sinusoid has been 
removed. This shows that the basic nonsinusoidal terms in the 
center-of-light velocity are well detected in our spectra. 

This solution has very large i, raising the interesting possi¬ 
bility of a pulsar eclipse. 


4. THE HIGH-ENERGY ORBITAL LIGHT CURVE 



4000 5000 6000 7000 

T c (K) 

Fig. 3. — Four-parameter ELC photometry/RV fit y 2 projected onto two 
planes; large dots indicate low y 2 - Square dots are for T c > 6000 K, con¬ 
sistent with the spectroscopic classification. Dotted boxes indicate locations 
of the various fits in Table|7](measured from finer grids). The "HiT” ELC-fit 
minimum consistent with the spectroscopy is shown in the inset box of the 
upper panel; here the point size scales inversely with y 2 , decreasing to 0 at 
Ay 2 /// = 3.0. Upper panel lines show the inferred pulsar mass. 
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Fig. 4.— Left: Top to bottom: gamma-ray pulse profile (with peak phase 
bins marked), orbital light curve for the peak phase bins (solid line) and full 
phase (dashed line, summed weights/8), and orbital 0.3-10 keV light curve, 
showing the broad minimum. The direct secondary eclipse phase is marked 
by the dashed lines, and error bars on the LAT-weighted count curves are 
computed as a = (1 + m?) 1 / 2 (see 2PC). Right: Likelihood ratio test 

of the weights in an orbital eclipse. The red line marks the measured TS, 
with 10 4 random realizations for d e = 40° (black histogram). Moderate 
significance is seen for eclipse widths 25°—45°. 


At our large i (“HiT”) fit minimum, we expect the com¬ 
panion to eclipse the pulsar for A 4>e = 25.4°. The radio¬ 
pulse emission is undetected for over half the orbit. However, 
the ionized companion winds of RB and black-widow pulsars 
generally prevent radio-pulse detection for sight lines passing 
well outside the companion Roche lobe, so this alone does not 
indicate true eclipse. 

We have re-examined the 17.0ks CXO ACIS expo¬ 
sure of PSR J2215+5135 (ObsID = 12466), starting on 
MJD 55697.1821 and covering 1.14 orbits. As noted by 
iGentile et al.l (12014ft . there is a clear shallow modulation with 
a minimum at pulsar superior conjunction 0 b ss 0.25 (Fig¬ 
ure 4). This shallow modulation might be best interpreted as 
variable viewing of an X-ray emitting intrabinary shock. 

A powerful intrabinary shock can also produce GeV LAT 
photons dXing & Wandl2015l) . However, since an apprecia¬ 
ble fraction of the object’s gamma-ray flux is pulsed, we 
can select the magnetospheric component and search for a 
true pulsar eclipse. Unfortunately, as for many RB MSPs, 
the wind obscuration and Pb fluctuations make radio tim¬ 
ing very difficult. At present, the best ephemeris available 
is that of lAbdo et al] OOBl 2PC), valid MJD 55346-55911, 
about a fifth of the LAT mission to date. We selected 10,951 
E = 0.1-30GeV Pass 7 reprocessed ‘Source class’ LAT 
photons within this date window and within 2° of the pul¬ 
sar. We compute weights Wi, the probability of being pulsar 
photons, with the LAT tool gtsrcprob using the lAcero et~afl 
(12015ft source spectrum parameters and the local background 
model; the mean source probability was (wi) = 0.037. The 
spin-phase 0s weighted light curve shows two peaks consis¬ 
tent with 2PC results (0s = 0.257 and 0s = 0.697 relative 
to the radio-pulse peak). If we extrapolate this ephemeris to 
the full LAT dataset, the peaks are lost. We select A 0 s = 0.1 
windows centered on the 2PC pulse peaks, to obtain events 
most dominated by magnetospheric emission. Other phases 
represent pulsed bridge flux, unpulsed flux from the intra¬ 
binary shock (if any), and unmodeled background emission. 
The middle panel shows the weighted binary light curve from 
the peak phase window; intriguingly, the minimum bin is at 
4>b = 0.25, pulsar superior conjunction. The full spin-phase 
light curve shows no strong eclipse. This suggests an addi¬ 
tional uneclipsed gamma-ray component. 


Binning clearly affects the minimum’s appearance, so we 
desire an unbinned test for the eclipse significance. For an 
eclipse width 6 e , we form the test statistic 

TS = Eilog [1 + Wi0 e /{ 1 - B e )] + E 2 log(l - w^, (2) 

where the first sum is over the uneclipsed window and the 
second is over the eclipse window. Exposure variations are 
already very small (~ 2% in the binned light curve), but 
the likelihood ratio test used here is insensitive to even these, 
since the variations are common to both the null and eclipse 
scenarios. To estimate the probability distribution of this 
statistic, we scrambled weights among the pulse peak pho¬ 
tons and recomputed TS 10,000 times. At 9 e = 25.4°, we 
find a chance probability of 0.080 (~ 1.4 a). However, the 
companion wind and photon opacity may widen the gamma- 
ray eclipse. Indeed, somewhat wider eclipses, to 9 e « 45°, 
show larger significance. The best detection is at 9 e « 40°, 
illustrated in the histogram at right, with a 0.7% (~ 2.45a, 
single trial) probability of a false-positive detection. This is 
intriguing, but hardly definitive. However, an improved ra¬ 
dio ephemeris should provide a factor of > 4 more photons 
and better pulse phase isolation, allowing a sensitive test for a 
magnetospheric (pulse phased) eclipse. 

5. CONCLUSIONS 

Our spectroscopic study of the RB MSP J2215+5135 does 
not support the inference, based on companion light-curve- 
shape fits, that it is a particularly massive pulsar. Instead, the 
results suggest a modest ~ 1.6M 0 NS mass. However, the 
present fit is based on a clearly inadequate heating model as¬ 
suming direct (X-ray) pulsar illumination. Of course, for all 
viable ELC model fits the heating flux is much larger than the 
Lx ~ 10 32 ergs -1 of the CXO observation. Indeed, as found 
for PSR J1311-3430 (Romani, Filippenko, & Cenko 2015), 
we find that the heating power into the companion solid angle 
is a large fraction of the spin-down flux (Figure 0 bottom); 
we argue that this indicates indirect heating by emission re¬ 
processed in an intrabinary shock. 

With the present direct heating model, the fits are inade¬ 
quate and the critical binary inclination i is not robustly de¬ 
termined. Accordingly, the fit masses are indicative only, and 
the statistical errors are not useful. There is a clear lesson 
here: with an inappropriate physical model, simple ELC fits to 
only the light-curve shape can deliver fit minima far from the 
true solution. Using the observed colors, as in ICARUS, cer¬ 
tainly helps, but spectroscopic constraints are also essential. 
Of course, amendment of the physical model should greatly 
improve the fit quality and may allow photometry-only solu¬ 
tions with high-quality data. 

Our temperature-constrained fits do suggest that the sys¬ 
tem is viewed nearly edge-on. If so, this presents the interest¬ 
ing possibility that the pulsar magnetosphere will be eclipsed 
by the companion. Limited by the radio ephemeris avail¬ 
able, present evidence for such eclipses is suggestive, but not 
definitive. However, additional radio data and further analy¬ 
sis of the Fermi LAT photons can produce a strong test of the 
eclipse hypothesis. If well measured, the width and spectrum 
of the pulsed flux eclipse can provide an important new tomo¬ 
graphic probe of the MSP magnetosphere and of the evapora¬ 
tive companion wind. 
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